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Abstract. We derive determinant expressions for the partition functions of 
spin-fc/2 vertex models on a finite square lattice with domain wall boundary 
conditions. 



0. Introduction 



In [1] , Korepin introduced the concept of domain wall (DW) boundary conditions 
for the six vertex (or spin-1/2) model on a finite lattice, and proposed recursion 
relations that fully determine the partition function in that case. In [2], Izergin 
obtained a determinant solution of Korepin's recursion relations. In this work, 
, inspired by Slavnov's inner product formula for higher spin models [3], we derive 

^ ' determinant expressions for the spin-fc/2 DW partition functions 1 , fc 6 N, using 

the fact that these models are related to the spin-1/2 model using fusion [4,5]. 

Basically, we show that fusion can be applied at the level of spin-1/2 partition 
functions to obtain spin-fc/2 partition functions, for any k G TV. More specifically, 
■ our result, in words, is that appropriate specializations of the rapidity variables in 

, Izergin's determinant expression for the spin-1/2 partition function on a kL x kL 

lattice (followed by suitable normalizations) yield determinant expressions for the 
^ c*| spin-fc/2 partition functions on an L x L lattice. 

In sections ^and [21 we briefly introduce the spin-fc/2 vertex models, and outline 
the fusion procedure. In|21and0] we motivate our result, then outline its proof. In 
El we give the spin-1 partition function as a specific example of our general result, 

S that has the benefit of allowing for an independent Izergin-type proof (which is not 
available for higher spin models) . In \§\ we derive the homogeneous limit of the 
, spin-fc/2 result, and comment on combinatorics of higher spin models. Finally, an 

appendix contains technical details. The presentation is elementary and (almost) 
self-contained. 



1. Vertex models 

Definitions related to vertex models. We work on a square lattice consisting 
of L horizontal lines (labelled from bottom to top), L vertical lines (labelled from 
left to right) and L 2 intersection points. 

We assign the z-th horizontal line an orientation from left to right, and a complex 
rapidity variable a^. We assign the j-th vertical line an orientation from bottom 
to top, and a complex rapidity variable yj. All rapidity variables are independent, 
unless specifically indicated to be otherwise. 

{x} is a set of rapidity variables {x\,--- ,xc}, where the cardinality of the set 
should be clear from context. 
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{x}j is a set {x\, x%, ■ ■ ■ , xc}, but with the Xi element missing. 

A fc-stack {x|fc} is a set of k variables of the form {x, x + 1, x + 2, • • ■ , a; + k — 1}. 

A bond is a line segment between two 
intersection points. 

A boundary, or extremal bond is a 

line segment, at the boundary of a line, 
attached to a single intersection point. 

In a spin-fc/2 models, we assign each 
bond k arrows, where k G {fc, k— 2, • • • , k 
mod 2}. All k arrows, on the same bond, 
point in the same direction. 

The k arrows on a bond define a spin 
variable on that bond. 
The magnitude of spin on a bond is 
k/2. 

The sign of spin on a bond is positive 
(negative) if the k arrows are oriented in 
the same (opposite) direction as (to) the 
rapidity flow in that bond. 

A vertex is the intersection point of the i-th horizontal line and the j'-th 
vertical line, together with the 4 bonds attached to it, and the arrows on them, 
a weight, toy, is a function assigned to a vertex, t>y, that depends on the difference 
of rapidity variables flowing through that vertex. In exactly solvable models, the 
weights satisfy the Yang Baxter equations [8]. 

Frequently used abbreviations. We frequently use the bracket notation [x] = 
sinh(Ax) (where A is a constant, 'crossing' parameter that characterizes the model) 
and the related product notation [x] TO = [a:] [a; — 1] • • • [x — m + 1], We also use the 
abbreviations x = x + 1 and iiy = — Xi + yj. 

fc/2xfe/2 

The partition function of a spin-fc/2 vertex model, on an L x L lattice, Z L xL , 
is a weighted sum over all configurations, that satisfy certain boundary conditions. 
The weight of a configuration is the product of the weights, toy, of the vertices t!y. 

(i) « fc/2 ({x},{y}) =e n «« 

vertices 

Conservation of spin flow. The vertex models, that we discuss in this work, 
conserve spin flow: In all vertices that have non-zero weight, the net incoming spin 
(measured with respect to rapidity inflow) equals the net outgoing spin (measured 
with respect to rapidity outflow). 

The spin- 1/2 model. There are six vertex types that conserve spin flow in the 
spin-1/2 model. They are shown in figure [21 below. For convenience, we label the 
vertices by their weights: an a vertex has weight a(x, y), and so forth. 

We do not need to distinguish vertices that share the same weight, except in one 
case (the c+ vertex) mentioned below. The six vertices of the spin-1/2 and their 
weights are shown in figure|21 The weights of every two vertices in the same column 
are equal and shown below them. 



x L ■ 



X2 
X\ 
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Figure 1. A finite square lattice, 
with oriented lines and rapidities. 
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a(x, y) = [u + 1] b(x, y) = [u] c(x, y) = [1] 

Figure 2. The six vertices of the spin-1/2 model and their weights, 
u = —x + y. 

The c+ vertex. The vertex with all 
arrows pointing inwards from left and 
right, and all arrows pointing outwards 
from above and below, plays a rather 
special role in this work. We refer to 
it as the c+ vertex. There is a unique 
c+ vertex in every spin-fc/2 model. The 
spin-2 c+ vertex is shown in figure [3] 

Domain wall (DW) boundary con- 
ditions. Consider the spin-1/2 model 
on a finite square lattice, and require 
that the boundary arrows have the same 
orientation as the arrow on the corre- 
sponding boundary of the c+ vertex: all 
arrows on the left and right boundaries 
point inwards, and all arrows on the 
upper and lower boundaries point out- 
wards. The c+ vertex is a DW configu- 
ration on a 1 x 1 lattice. 



The Izergin-Korepin (IK) determinant expression for the spin-1/2 partition 
function [1,2] is 

L 

11 [~Xi + y 3 + l]a 

( 2 , (m.m) - + «,H-„ K 2 "' 2 ) 

l<i<j<L 

where 
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Definitions related to partition functions, matrices and normalizations. 

k /2xk /2 

Z L xL is the partition function of the spin-fc/2 model on an L x L lattice. 

M^ 2 £ k ^ 2 is the L x L IK matrix, but with a k x k block structure, as will be 
explained below. 

^k^Lxkl% IS the kxk ij-th block of M^ 2 ^^ 2 . It depends on the rapidities {xi, yj}. 
B k * k is the ij-th kx k block of the lattice. 

^kLxkL 2 * s tne normalization function of M k ^ k ^ 2 that sets the weights of the 
spin-fc/2 vertices to [k)k- 

Spin-fc/2 models. We will not need the weights of spin-fc/2 model in all generality. 
They can be deduced from those of the fused elliptic height models, [5], as follows. 

Spin-fc/2 vertex weights from elliptic height weights. 

§1. Take the trigonometric limit: set the elliptic nome, which appears in the 
weights of elliptic models, q — > 0. This reduces the weights to ratios of prod- 
ucts of trigonometric functions (for pure imaginary values of the crossing 
parameter). 

§2. Take the vertex limit: set the height shift parameter, that appears in the 
weights of height models, ( — ► ±oo. This eliminates dependence on the 
height variables. 

§3. Symmetrize the resulting weights so that the weights of the c-type vertices 
are equal. 

2. Fusion 

As the spin-fc/2 vertices are obtained from the spin- 1/2 vertices using fusion, we 
wish to recall how fusion works, following [5] 2 . 

Definitions related to boundaries. A boundary of length L is a set of L par- 
allel extremal bonds. A vertical boundary consists of horizontal extremal bonds. A 
horizontal boundary consists of vertical extremal bonds. Notice that, in our defi- 
nition, a boundary cannot be closed. A closed boundary will consist of a sequence 
of horizontal and vertical boundaries. 

A cr-configuration is an arrangement of spins on a boundary, with total spin a. A 
boundary of length L, in a spin-fc/2 model, has a 6 {fcL/2, fcL/2 — 1, • • • , — kL/2}. 
A CT-set is the set of all ^-configurations on a boundary. 

A cr-representative of a <r-set, is any single uniquely defined configuration in that 
set, that we select to represent the entire set. In this work, we choose that to be 
the configuration with spins ordered in terms of their values, with larger (more 
positive) spins to the left of (lower than) smaller (less positive) spins in the case of 
horizontal (vertical) boundaries. 



An inflow (outflow) boundary is one 

that rapidity variables flow into (out of), 
as seen from the inside of the region that 
it bounds. 



6/2 



a/2 



7 /2 



0/2 



Figure 5. Spin values of a gen- 
eral vertex in a spin-k/2 model. 



Although [5] studies fusion in the context of elliptic height models, it is convenient to start 
from there, as their exposition is explicit. 
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Fusion procedure. Following [5], to compute the weight of the generic spin-fc/2 
vertex shown in figure^] where (a, (3, 7, 5) £ {k, k— 1, • • • , —k} and a+f3 = j+S, we 
start from the set of all spin-1/2 configurations on a k x k lattice, with boundaries 
that match in total spin values those of the vertex that we wish to produce (the 
right boundary has total spin a/2, etc), and proceed as follows. 

§1. Set {xi, x 2 , • ■ ■ , x k } to {xi|fc}, and set {yi, y 2 , ■ ■ ■ , 2/fc} to {yi|fe}. x x and 
?/i will be the rapidities of the resulting spin-fc/2 vertex. 

§2. Sum over all a/2-configurations on the left (inflow) boundary. 

§3. Sum over all /3/2-configurations on the lower (inflow) boundary. 

§4. Take the 7/2-configuration on the right (outflow) boundary to be the unique 
7/2-representative. No summation over configurations is performed. 

§5. Take the 5/2-configuration on the upper (outflow) boundary to be the 
unique (5/2-representative. No summation over configurations is performed. 

§6. Normalize the result 3 , by dividing with 



fc-i 

(3) N^l xk/2 (x uyi ) = Y[[-x 1+yi + p} k ^ 

Remarks. To perform fusion, in our convention, inflow boundaries are summed 
over, while outflow boundaries are set to representative configurations. Further, we 
could have set the outflow boundaries to any cr-configuration with the correct net 
spin <7. However, using the Yang-Baxter equations, one can show that the result is 
independent of the choice [5] . 

3. SPIN-fc/2 PARTITION FUNCTIONS: MOTIVATION 



Suppose we wish to obtain the spin-2 
c+ vertex of figure [21 Following the fu- 
sion procedure, we need to consider the 
4x4 DW spin-1/2 partition function, 
shown in figure |21 

But this case is very simple: Because 
of the boundary conditions, all cr-sets 
have exactly one configuration each, and 
we do not need to sum over left and lower 
configurations, or choose any right and 
upper ones. All we need to do is to set 
the rapidities to the right values, and 
normalize suitably. 



X\ 



X3 



X 2 



X\ 



Vl 2/2 2/3 2/4 

Figure 6. 4 x 4 spin-1/2 parti- 
tion function. 



Now, suppose we do not wish to fuse 
all the way down to the c+ vertex, which 
is, the lxl DW spin-2 partition func- 
tion, but only half way to the 2x2 DW 
spin-1 partition function. 



•1'2 



2/1 2/2 

Figure 7. 2 x 2 Spin-1 partition 
function. 



3 Our normalization is different from that of [5]. We choose to normalize the c+ vertex of 
the spin-fc/2 model to [k]^. In [5], the c+ vertices are normalized to 1 (up to phases), in the 
trigonometric vertex limit that we are interested in. 



6 



CARADOC, FODA, AND KITANINE 



Partial fusion. It seems plausible that all we need to do in this case is to 2-stack 
the rapidity variables, and normalize suitably. In other words, we need partial 
fusion as follows. 

§1. Consider the spin- 1/2 model on a kL x kL lattice with DW boundary 
conditions. 

§2. Set {xi, X2, ■■■ , xul} to L fc-stacks {{xi|fc}, {x 2 |/c}, ••• , {xl|A;}} and 
{yx, V2, ■■■ , Vkh} to L /c-stacks {{yi|fc}, {yz\k}, • • • , {y L |fc}}. Under this 
restriction of variables, the Ml^u/ 2 IK matrix, that we started with, is 
now denoted by Mj,0*^ 2 . 

§3. Normalize so that the weight of the spin-fc/2 c+ vertex is a constant 4 by 
dividing with 

<r /2 ({x},m) = n n l H*+w+p]*-x 

l<i,j<L p=0 

Following the above procedure, we obtain the following expression for the spin- 
fc/2 partition function 



(4) z k L / x T /2 ({x},{y}) = 

k 

n n \~ x i +yj +p\k+i 

l<i,j<L p=l 

k-1 k-1 

n n [~ x i + x j + p\k n + v* + ?]* 

l<i<j<£p=0 p=0 



Equation 0| is our main result. In the next section, we show that partial fusion, 
as outlined above, works, and leads to spin-fc/2 partition functions, with no missing 
or unwanted configurations. Technical details of how equation 0] is obtained are in 
the appendix. 

4. SPIN-fc/2 PARTITION FUNCTIONS: PROOF 

We wish to show that, starting from a weighted sum over spin-1/2 configurations, 
on a kL x kL lattice, and dividing the lattice into L 2 k x k blocks, we can fuse 
these blocks one by one, and obtain a weighted sum over spin-fc/2 configurations 
with the correct spin-fc/2 weights. 

In particular, we also wish to show that This partial fusion procedure leads to 
all required spin-/c/2 configurations, and no more. In fact, it will turn out that this 
procedure is bijective in the sense that every step is reversible. 

The following is an outline of the proof, together with a simple running example. 

4.1. Outline of proof. 

§1. Consider a DW spin-1/2 model on a kL x kL lattice. As an example, we 
take k = 2 and L = 3. 



The rationale of normalization is to put the result in a practical, recognizable form, and in 
particular to avoid that the weights have common, spurious poles or zeros. 
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§2. Set the horizontal rapidities into L fc-stacks of the form {{xi|fc}, {x 2 |fc}, 
• • • {xl|/c}}, and similarly for the vertical rapidities. In our example, we 
obtain {{zi|2},{a;2|2},{:c3|2}} and {{xx|2}, {x 2 |2}, {x 3 |2}} 

§3. Divide the lattice into L 2 k x k blocks, B k * k , where i is the block row 
index, j is the block column index, and {i,j} G {1, 2, • • • , L}. Each block 
has one independent horizontal rapidity Xi, and one independent vertical 
rapidity yj. In our example, we obtain nine 2x2 blocks. 

§4. Order the blocks, firstly in terms of row position: blocks with a smaller i 
precede those with a larger i, then in terms of column position: for equal i 
indices, blocks with a smaller j precede those with larger j. In our example, 

+ 1^ ™-rW ;= / R 2x2 R 2x2 R 2x2 H 2x2 H 2x2 R 2x2 R 2x2 R 2x2 R 2x2 X 

the order is { B lx , B 12 , B 13 , B 21 , B 22 , B 23 , B 31 , B 32 , B 33 } 



X3 
X2 

x 2 



X\ 
X\ 



yi yi 2/2 2/2 2/3 2/3 

Figure 8. Spin-1/2 partition function with 2- 
stacked rapidity variables, z = z + 1 . 

l/2xl/2 

§5. Consider the partition function Z k ' Lxk £ as a sum over products of two 
partition functions: that of B^ k , and that of the rest of the lattice, i?n, 
that is 

(5) E < fe *n 

{hi,vi} 

where hi and vi stand for the horizontal and vertical common spin bound- 
aries. The sum is over all configurations on the common boundaries of 
B h * k and i?n. 

§6. Each B k ^ k term in the above sum is DW on an inflow boundary, and has 
a single cr-configuration on an outflow boundary. Since a DW boundary 
condition corresponds to a cr-set with a single cr-configuration, the inflow 
boundaries are (trivially) summed, while the outflow boundaries are fixed 
to a certain cr-configuration. But, from fusion, all such partition functions 
are equal to the one with cr-representatives on the outflow boundaries. This 
allows us to simplify the sum in the previous equation to 

(6) « 1/2 = £(< fe £*n) 

where the first (right most) sum is over all cr-configurations, in an allowed 
cr-set on the inflow boundaries of Ru, B k f k has cr-representatives on the 
outflow boundaries, and the second (left most) sum is over all cr-sets on the 
common boundaries. Notice that we need to be clear about what is being 
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summed over in each sum. This is because, before fusion, all elements in all 
allowed d-sets are summed over, while, after fusion, only cr-representatives 
are summed over. 



x 3 
xz 

X2 

x 2 



hi 



yi 



hi 



Vl 



2/1 



2/2 2/2 2/3 2/3 



Figure 9. Detaching a 2 x 2 spin-1/2 block, z = z + 1. 

§7. Using fusion to write B^ k as a spin-&;/2 vertex, we end up with a sum 
over products of a 1 x 1 spin-/c/2 non-DW partition function and a spin-1/2 
non-DW partition function (the original lattice minus B^* ). 



(7) 



7 l/2xl/2 
J kLxkL 



k/2 7 l/2xl/2 
/ , ^{kLxkL)-( 



ixl)) 



^3 
£3 
X 2 
X 2 



Xi- 



2/1 



2/2 2/2 



2/3 2/3 



Figure 10. Fusing a detached 2x2 spin-1/2 Wocfc into a 1 x 1 
spin-1 partition function. 2 = 2 + 1. 



Next, we consider the next ranking block in Z, kLxkL ,_, lxl , 
and write 



1/2x1/2 namely i?***, 



(8) 



7-1/2x1/2 
7 0LxA:L)-(lxl) 



Enfexfe 7 l/2xl/2 
-°12 ^(kLxkL)-(2xl) 
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(9) 



7 l/2xl/2 
J (kLxkL)-(lxl) 



E 



1/2x1/2 
(kLxkL)- 



(2x 



Using fusion to re- write B^ k as 
a spin-fc/2 vertex, combining the 
above results, and summing over the 
common boundaries of the two 1 x 
1 spin-fc/2 partition functions, we 
can write the initial spin- 1/2 non- 
DW partition function as a sum over 
products of two objects: a spin-fc/2 
non-DW partition function consist- 
ing of 2 vertices, and the remaining 
spin- 1/2 lattice. 

It should be clear from the above 
that we are fusing the initial spin- 
1/2 lattice, one block at a time, to a 
spin-fc/2 lattice. It should also be 
clear that this can be done block 
by block, that the procedure is re- 
versible, and that we end up with the 
desired L x L DW spin-fc/2 model. 

The DW boundary conditions of 
the final, spin-fc/2 configurations fol- 
low from the DW boundary condi- 
tions of the initial spin-1/2 configu- 
rations. 



£3 



X-2 



■f'2 



Xl. 



2/3 V3 



V2 



m 

Figure 11. Detaching, fusing the second 
2x2 spin-1/2 block, fusing it to form a 
spin-1 vertex, then attaching the latter to 
the first spin-1 vertex, z = z + 1. 



This concludes our outline of the proof of equation^] from which the reader can 
recover a formal proof if necessary. 



5. Example: Spin-1 model 

Consider the spin-1 (19-vertex, or Zamolodchikov-Fateev model [9]) model con- 
structed by fusion 5 . The vertex weights, which can be computed explicitly using 
fusion, are listed (for example) in [10]. In the following, we list only those vertices 
that we need and their the weights. Vertices that are related to those listed, by 
inverting of all arrows, have the same weights. 



^There are many 19-vertex models. The Zamolodchikov-Fateev model that we are considering 
here is only on of these. Two others are listed in [10]. 
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B=[u-l][u] Y=[u][2] E = [u][u+1] 

Figure 12. A subset of the vertices of the spin-1 model and their 
weights, u — —x + y. 
Re- writing equation^] for k = 2, in terms of vertex weights of the spin-1 model, 
we obtain the following expression for the spin-1 partition function 



(10) Z}£ L ({x},{y}) = 

II A(-Xi + yj)E(-Xi + yj)B(- Xi + yj) 

l<i,j<L 

11 E (~ x i + Xj)B(-Xi + Xj)E(yi - y J )B(y l - yj) 

xdet {M^ 2L ) 



= W-Xi + Vj ), M|/_ 2 li2j - 1/B(- Xi + Vj ), 

(11) 

Mlf 2l _! = l/A(- Xi + yj ), M^ 2 2l = 1/E(-Xi + y s ) 

Independent check. The determinant expression of the spin-1 partition function, 
obtained above, allows for an independent check of our general result, in the sense 
that one can take the determinant expression as a conjecture, and show that it is 
correct. 

Just as in the spin-1/2 case [1,2], one needs to show that the LHS of equation 
I1UI (the partition function) has certain properties that uniquely determine it com- 
pletely, then show that the RHS (the proposed determinant expression) satisfies 
the same properties. An outline of the basic steps is as follows. 

The LHS of equation 1101 The following properties fully characterize the parti- 
tion function 

Symmetry Using the Yang Baxter equations, one can show that the partition 
function is a symmetric function in {x} and in {y}. 

Degree One can easily show that, on any extremal row or column, there is exactly 
one rapidity independent vertex (namely a c+ vertex), or two vertices that are 
trigonometric polynomials of degree 1 each, while all other vertices are of degree 
2. This means that the partition function is a trigonometric polynomial of degree 
2L — 2 in the rapidity variable in that row or column, and by symmetry in every 
other variable. 

Recursion relations Below, we will show that the partition function satisfies 2L 
recursion relations in each rapidity variable, which is more than we actually need. 
The initial condition By construction, the lxl DW partition is a c+ vertex. 
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The RHS of equation 1101 Next we show that the RHS satisfies the same prop- 
erties as the LHS. 

Symmetry By direct calculation, one can show that the RHS is a symmetric 
function in {x} and in {y}. 

Degree Naive power counting shows that the RHS is a degree 2L trigonometric 
polynomial in any rapidity variable x. However, taking the limit x — > oo, for real 
crossing parameter, one can show explicitly that degree 2L terms cancel, while 
there are no degree 2L — 1 terms, so that the RHS is a degree 2L — 2 trigonometric 
polynomial. 

Recursion relations The RHS satisfies the same recursion relations as LHS, as 
will be shown below. 

Initial condition By direct calculation, the RHS for L — 1 reduces to the weight 
of the C+ vertex. 



Recursions from upper left corner. Due to the boundary conditions, the only 
vertices that are allowed at the upper left corner are A(—x\ + yi), X{—x\ + y\) 
and C(—x\ +yx). 

Setting —xi + y\ + 1 — 0, we obtain A(— 1) = X(— 1) = 0, which freezes all 
vertices on the top row and first column and leads to the recursion relation 



(12) Z)£ L [{*},{y}\x 1 =y 1 + l) 
[1][2] 



r l 



Y[ B(-xi + Vj-xiBi-Xj + yi) 
U= 2 



^-l)x(L-l) (Ml,Ml) 



Given the symmetry in vertical rapidities, we get the same relations for Xi — 
yj + 1, for all i and all j, so we have L recursion relations for each rapidity variable. 



Recursions from upper right corner. Due to the boundary conditions, the only 
vertices that are allowed at the upper right corner are B{—X\ + yi,), Y{—x\ + yi) 
and C(-xi + y L )- 

Setting — x\ + yi = 0, we obtain B(0) = Y(0) = 0, and only C(0) survives at 
the corner, freezing all the vertices on the top row and last column. The remaining 
[L — 1) x {L — 1) lattice has once again DW boundary conditions, and we obtain 
the recursion relation 



(13) Z\\\ ({x},{y} 
[1][2] 



%i = Vl 



Y[A(-xi + yj )A(-*i 

U=2 



VL) 



7 lxl 

J (L-l)x(L- 



-l) ({x}i,Ml) 



Given the symmetry in vertical rapidities, a similar relation, Xi = yj, can be 
written for all i and j, and we have L recursion relations for each rapidity variable. 
Thus we have altogether 2L recursion relations for each variable, which are sufficient 
to completely determine the partition function as a trigonometric polynomial of 
degree 2L — 2 in that variable, just as Izergin's proof. 

As the RHS of equation I1UI satisfies the 2L required recursion relations and is a 
polynomial of degree 2L — 2, in every rapidity variable, and given that the initial 
condition is satisfied, we conclude that it coincides with the LHS of equation 1101 
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Comments. The above Izergin type proof does not extend to higher spin models, 
beyond spin-1. The reason is that, for k > 2, the degree of the polynomials that 
we need to determine is higher than the number of available recursion relations 6 . 



6. The homogeneous Limit 
Taking the homogeneous limit of the spin-fc/2 partition function, following the 



footsteps of [2], is straightforward. For convenience, we re- write M ; 



k/2xk/2 
kLxkL,(i,j) 



as 



M, 



k/2xk/2 _ 
kLxkL,(i,j) 

' <t>{-Xi +Vj) 
<t>{-Xi + y.j - 1) 



, (f)(-Xi + yj - k + 1) 



{-Xi + y 3 + 1) 
i-Xi+yj) 



+ yj + k-l) ) 



<f>(-Xi + Vj) 



where 4>{x) = l/(\x][x + 1]). Let x\ = x and consider the limit X2 — > x. The first 
block row M^tJiIt ?j remains unchanged apart from replacing x\ with x, while 



M; 



kL x kL,lj 

k/2xk/2 , 

kLxkL,2j becomes 



M, 



k/2xk/2 _ 
kLxkL,(2,j) ~ 

4>(-X 2 +Vj) 
<t>(~X2 + Vj - 1) 



-X2 + Vj — k + 1) 



(-X2 + Vj + 1) 
{-X2 +Vj) 



<t>(-X2 + Vj + k ~ 1) ^ 



-^2 + Vj) 



Taylor expanding each term as x^ — > x, to first order, and subtracting the first 
block row from the second 



<j)'{-x + yj) <j>' {-x + + 1) ... ^(-x + yj + k-1) ' 

tfi-x + yj-l) 4>'{-x + yj ) : 

^(-jB + y^-jfe + l) ... + 

where <j^ n > (x) represents the n th derivative of 4> with respect to its argument. Note 
that the [— x^ +x] k term in the denominator of equation 01 cancels exactly with the 
overall factor (x2 — x) k as X2 — > x. 

Taylor expanding as Xj — » a;, up to n-th order, and successively eliminating 
terms by subtracting multiples of previous block rows and taking out common 
factors, M^xkL% becomes 



(x 2 - x) k 



This assumes that we are willing to consider only Lagrange interpolation in determining the 
polynomials under consideration. One can consider more elaborate interpolations, such as Hermite 
interpolation, but then things become very complicated. 
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(xi - x)W-V ((i-l)l) 
^-^(-x + yj-l) 



^(-x + yj + l) 



»-i)r- 



(-x + yj + k-l) 



^-^(-x + yj-k + l) ... fi^-V (-x + yj) 

As before, the denominator in equation ^contributes a factor of [— £j + a;] fc ^~ 1 ^ ) 

which cancels with the above coefficient as Xi — > x, leaving a factor of f (i — 1)! J , 
and semi-homogeneous spin-fc/2 L x L partition function becomes 



(14) 



L k 

n Yi[-x+yj+v]k 



fc-l , V 

n (-M 2 ) 

i=i v y 



(fc-i)(L 2 -L)/2 L-l 



n [*!] n n 

=1 v J l<i<j<Lp=l 



[-Vi + Vj +P]k 



With M, 



k/2xk/2 
kL x kL,ij 



given by the above matrix. Eauation ll4l is the partition function 
for a lattice where the horizontal rapidities are homogeneous but the vertical ra- 
pidities are still distinct. Applying similar arguments to the vertical rapidities, and 
combining results, we obtain the following expression for the homogeneous L x L 
spin-/c/2 partition function 



(15) 



2 



fc/2xfc/2 
LxL 



where M%?£*! 2 is 



II [-x + y +p]k 



det {Mt' L \\T) 



k-i f x 

n \-w) 



(k-i){L 2 -L) L-l 



n («) 

=i v y 



2/, 



0(*+i-2)( u ) 
0(*+i-2)( u _ i) 



(i+i-2)( u + i) 
(«'+i-2)( u ) 



(^- 2 )(u + fc-l) ) 



(i+i-2) 



(«) 



Remarks on combinatorics. Do the higher spin determinants lead to interesting 
combinatorics, as in the spin-1/2 case? It is not difficult to show that there is 
a simple bijection between spin-fc/2 DW configurations and extended alternating 
sign matrices ASM's with entries in {0, ±1,- • ■ ,±fc}, and conditions that naturally 
extend those of the usual ASM's [11]. 

However, one can also easily check that, unlike the spin-1/2 case, for k > 2 
there is no choice of crossing parameter and rapidity variables such that all weights 
are equal (even up to phases). This rules out 1-enumerations (but not weighted 
enumerations). This conclusion is substantiated by direct numerical enumerations 
(for k = 2), that lead to numbers that cannot be expressed as simple products [11]. 
We hope to return to these issues in a separate publication. 
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Appendix: Technical details 
We start from the partition function of the spin- 1/2 model on a kL x kL lattice 



(16) 



7 l/2xl/2 
J kLxkL 



EI [~ x i + V] + 1] [~ x i + Vi\ 
l<i,j<kL 



n 



l<i<j<kL 



n i-vi + Vj\ 



det 



( M l/2X1/2^| 



l<j<i<kL 



To obtain the fused partition function, we proceed in two steps. Firstly, we set 
the rapidities to suitable values, then we normalize the result. Setting the rapidities 
to {{xi|/c}, {x 2 |fc}, ••• , {x L |/c}}, and {{yi|fe}, {y 2 |fc}, ■•• , {yL|fc}}, and using 
Uij = —Xi + yj, and Vij = —xi + Xj, we obtain 



(17) n [uij] -» 

l<i,j<kL 

Yl [uij -k+l][uij -k + 2] 

l<i,j<L 



2 [ Uij ] k ■ ■ ■ [uij +k- 2] 2 [My + k - 1] 



(18) [] fal 

l<j<i<fe_L 

JJ — fe + 1] [wy - fc + 2] 2 • ■ ■ [«<,•]* ■ • • [vij +k — 2} 2 [ VlJ +k-l] 

l<j<i<L 

p=l 

where the last factor in eciuation ll8l comes from diagonal terms such as [— a^+xi] — > 
[— (xi + l) + a;i] = [— 1]. Using eauationslT7land [T8l the RHS of eauationllBlbecomes 



(19) 



n 

l<i,j<L 



' k 

n [-Xi+vj+p\k+-i 



fk-i 

I! [~ x i +Vo +P\k-1 
{p=o 



det (M^f) 



(k-l 




2L 








fe-i 


n 


i-p]p 




n Til-*- 


hxj +p] k 




n n [-vi + v 3 










> 




t l<i<j<L p=0 



To normalize, we divide by a factor so that the weight of the c+ vertex in the 
spin-fc/2 model is [k]k, and multiply by a factor that accounts for the change in 
symmetry structure of the determinant. Together, these are 



iri-dp 



2L 



fc-1 



n n [- x i+yj +p\k-i 

K l<i,j<L p=0 

and we end up with the general result in equation 01 
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